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m ; Abstract 

© ; We develop a simple model for a self-gravitating spherically symmetric relativistic star which 

begins to collapse from an initially static configuration by dissipating energy in the form of radial 
Q-c 

heat flow. We utilize the model to show how local anisotropy effects the collapse rate and thermal 
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p,i behaviour of gravitationally evolving systems. 
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I. INTRODUCTION 

In cosmology and astrophysics, there exist many outstanding issues relating to a dynam- 
ical system collapsing under the influence of its own gravity. In view of Cosmic Censorship 
Conjecture, the general relativistic prediction is that such a collapse must terminate into 
a space-time singularity covered under its event horizon though there are several counter 
examples where it has been shown that a naked singularity is more likely to be formed (see 



w 



and references therein). In astrophysics 



long been very much speculative in nature 



;he end stage of a massive collapsing star has 



2] . From classical gravity perspective, to get a 



proper understanding of the nature of collapse and physical behaviour of a collapsing sys- 
tem, construction of a realistic model of the collapsing system is necessary. This, however, 
turns out to be a difficult task because of the highly non-linear nature of the governing field 
equations. To reduce the complexity, various simplifying methods are often adopted and 
the pioneering work of Oppenheimer and Snyder [31] was a first step in this direction when 
collapse of a highly idealized spherically symmetric dust cloud was studied. Since then, 
various attempts have been made to develop realistic models of gravitationally collapsing 
systems to understand the nature and properties of collapsing objects. It got a tremendous 
impetus when Vaidya[4] presented a solution describing the exterior gravitational field of a 
stellar body with outgoing radiation and Santos[5| formulated the junction conditions join- 
ing the interior space-time of the collapsing object to the Vaidya exterior metric [4J]. These 
developments have enabled many investigators to construct realistic models of gravitation- 
ally evolving systems and also to analyze critically relevance of various factors such as shear, 
density inhomogeneity, local anisotropy, electromagnetic field, viscosity etc., on the physical 
behaviour of collapsing bodies|6l-ll9l I2ll-l4ll.l44j-l46|. In the absence of any established theory 



governing gravitational collapse, such investigations have been found to be very useful to 
get a proper understanding about systems undergoing gravitational collapse. 

The aim of the present work is to develop a simple model of a collapsing star and in- 
vestigate the impact of pressure anisotropy on the overall behaviour of the collapsing body. 
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42|, 



Anisotropic stresses may occur in astrophysical objects for various reasons w 

phase transition, density inhomogeneity, shear, electromagnetic field etc. J9|, |42j, |47| 

it has been shown that influences of shear, electromagnetic field etc. on self-bound systems 

can be absorbed if the system is considered to be anisotropic, in general. Local anisotropy 



has been a well motivated factor in the studies of astrophysical objects and its role on the 
gross features of static stellar configurations have been investigated by many authors (see 
for example, Ref. [9l. llOl l43l. I48l-l5l| and references therein). For dynamical systems, though 
pressure anisotropy is, in general, incorporated in the construction, very few works have been 
reported till date where impacts of anisotropy have been discussed explicitly |llHl3l |20| . |23 |. 



Appropriate junction conditions for an anisotropic fluid collapsing on the background space- 
time described by the Vaidya metric has been obtained in 24|. Considering a spheroidal 
geometry of Vaidya and Tikekar 54| , Sarwe and Tikekar 4J] have analyzed the impact of ge- 
ometry vis-a-vis matter composition on the collapse of stellar bodies which begin to collapse 
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46| 



from initial static configurations possessing equal compactness. Sharma and Tikekar 
have investigated the evolution of non-adiabatic collapse of a shear-free spherically symmet- 
ric object with anisotropic stresses on the background of space-time obtained by introducing 
an inhomogeneous perturbation in the Robertson- Walker spacetime. 

In the present work, we have developed a model describing a shear-free spherically sym- 
metric fluid distribution radiating away its energy in the form of radial heat flux. The star 
begins its collapse from an initially static configuration whose energy-momentum tensor 
describing the material composition has been assumed to be anisotropic, in general. To 
develop the model of the initial static star, we have utilized the Finch and SkeaJ52j ansatz 
which has earlier been found to be useful to develop physically acceptable models capable of 



describing realistic stars in equilibrium 
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57] . The back ground space-time for a static 



configuration for the given ansatz has a clear geometrical interpretation as may be found in 
Ref. 56] . By assuming a particular form of the anisotropic parameter, we have solved the 
relevant field equations and constructed a model for the initial static stellar configuration 
which could either be isotropic or anisotropic in nature. The solution provided by Finch 
and Skea|52|] is a sub-class of the solution provided here. Since, the solution presented here 
provides a wider range of values of the anisotropic parameter, it enables us to examine the 
impact of anisotropic stresses on the evolution of a large class of initial static configurations. 

Our paper has been organized as follows: In Section [Til we have presented the basic 
equations governing the system undergoing non-adiabatic radiative collapse. In Section [TTTl 
by assuming a particular anisotropic profile, we have solved the relevant field equations 
to develop a model for the initial static star. In Section IIV} by stipulating the boundary 
conditions across the surface separating the stellar configuration from the VaidyaJ4| space- 



time, we have solved the surface equation which governs the evolution of the initial static 
star that begins to collapse when the equilibrium is lost. In Section |V] we have analyzed 
the impact of anisotropy on gravitationally collapsing systems by considering evolution of 
initial static stars which could either be isotropic or anisotropic. Impact of anisotropy on the 
evolution of temperature has been analyzed in Section IVI1 Finally, some concluding remarks 
have been made in Section IVIII 

II. EQUATIONS GOVERNING THE COLLAPSING SYSTEM 



We write the line element describing the interior space-time of a spherically symmetric 
star collapsing under the influence of self-gravity (in standard coordinates x° = t, x 1 = r, 
0, x 3 = 



x 



as 
ds 2 



-A 2 Q (r)dt 2 + f 2 (t)[B 2 (r)dr 2 + r 2 (d6 2 + sin 2 6d(f)% (1) 

where, A (r), B (r) and /(£) are yet to be determined. Note that in Eq. (JT]), if we set 
f(t) = 1, the metric corresponds to a static spherically symmetric configuration. 

We assume that the matter distribution of the collapsing object is an imperfect fluid 
described by an energy-momentum tensor of the form 



T a p = (p + Pt)u a up + p t g a p + (p r - Pt)XaXf} + QaUp + qpu c 



(2) 



where, p represents the energy-density, p r and p t , respectively denote fluid pressures along 
the radial and transverse directions, u a is the 4-velocity of the fluid, x a is a un it space-like 
4-vector along the radial direction and q a = (0, ^,0, 0) is the heat flux vector which is 
orthogonal to the velocity vector so that u a u a = — 1 and u a q a = 0. 

The Einstein's field equations governing the evolution of the system are then obtained as 
(we set G = c = 1) 
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In Eqs. ©-(JBj), a 'prime' denotes differentiation with respect to r and a 'dot' denotes dif- 
ferentiation with respect to t. Making use of Eqs. (Jlj) and (jSJ), we define the anisotropic 
parameter of the collapsing object as 



A(r,t) = 87r(p r - p t ) = — 



Aq ^o | -^o | A B' 1 
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We rewrite Eqs. (j3j) 
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where p s , (p r ) s an d {pt)s denote the energy-density, radial pressure and tangential pressure, 
respectively of the initial static star. 

III. INTERIOR SPACE-TIME OF THE INITIALLY STATIC CONFIGURATION 



In our construction, we assume that an initially static star (with /(£) = 1 in Eq. ([[])), 
described by metric potentials Aq(t), Bq(t) and anisotropy A s (r), starts collapsing if, for 
some reasons, it loses its equilibrium. To develop a model of the initially static configuration, 
we first assume that the anisotropic parameter is separable in its variables so that A(r, t) = 
A s (r)/p(t). Eq. © then reduces to 
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which is independent of t. Eq. (11 ip can only be solved if any two of the unknown functions 
(A , Bq and A s (r)) are specified. To develop a physically reasonable model of the initial 
static configuration, we first utilize the Finch and Skea 52| ansatz for the metric potential 
Bq given by 

/ l + i\ (12) 



B 2 (r) 



B 2 ) 



where R is the curvature parameter describing the geometry of the configuration. In the 
static case, the ansatz (fT2|) has a clear geometric characterization and has been found to 



yield realistic models for compact stellar objects[55N57]. Substituting Eq. (TT2|) in Eq. (fiTj) . 
we obtain, 



x- 



cPAo _ 2 n-x 2 \ dA 

dx 2 \ R 2 x 5 ) dx 

where we have used the following transformation: 
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A, = 0, 
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1 + R 2 )- 



(14) 



To solve Eq. (1131) . we assume a particular anisotropic profile 

a(x 2 - 1)(2 - x 2 ) 



A s (x) 



R 2 x % 



(15) 



where a is the measure of anisotropy. The motivation for choosing the particular form of 
the anisotropy parameter are the following: (1) It is physically reasonable as the anisotropy 
vanishes at the centre (r = 0, i.e., x — 1) as expected and (2) it provides a solution of 
Eq. (fTBl in closed form. Note that a = corresponds to an initial static star which is 
isotropic in nature. Substituting Eq. ( fl5l) in Eq. ( fl3l) . we get 



d 2 Ar 



2dA 



dx 2 x dx 



+ 



a(x 2 -2) 



;ir 



A n = 0, 



(16) 



whose solution is obtained as 



Mx) = Px 3/2 Ji 



v / 9 = 8a 1 ' 
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V-l + a) + Qx 3/2 Yi V g =M (-ixV-l + a) 



(17) 



where P and Q are integration constants; Ji v ^psa(— ix\J — 1 + a) is the Bessel function of 
first kind of order |\/9 — 8a and yi^ 9 _ gQ , (— %X\J— 1 + a) is the Bessel function of second 



kind of order |\/9 — 8a. It is obvious that solution is valid for a < 1. At a = 1, the 
Bessel function encounters a singularity and, therefore, we shall deal with the a = 1 case 
separately. 



Special cases: 

• Case I : a = (A,(r) 
Eq. (TT31) reduces to 



0), i.e., initial static configuration is isotropic in nature. 



d 2 A 2 dA 



+ A = 0, 



(ia; 2 x <ir 
6 



(18) 



whose solution is found to be 

A (x) = [(C - Dx) cosx + (Cx + D) sinx], (19) 

where, C and D are integration constants. This particular solution has been found 



previously in Ref. 52J. However, the solution ffT9|) can be obtained directly from the 
solution (TT71) by substituting a = and setting C = — y~Q, and L> = w |P. 

Case II : a = 1 (A s (r) ^ 0). 

Eq. (fl6l) reduces to 

£• _ 2£> + ^ = , (20) 

dx A x ax x A 



whose solution is given by 



A (x) = Ax + Bx 2 , (21) 



where A and B are integration constants 55[]. Unfortunately, the above solution can 



not be regained from the general solution (TT7|) due to the properties of Bessel functions 
and should be treated separately 

We thus have a model for an initially static stellar configuration which could either be 
isotropic or anisotropic in nature. After loss of equilibrium, the initially static star starts 
collapsing and to generate a solution of the subsequent evolving system, we need to determine 
f(t). This will be taken up in the following section. 

IV. EXTERIOR SPACE-TIME AND JUNCTION CONDITIONS 

In our construction, evolution of the collapsing object is governed by the function f(t) 
which can be determined from the boundary conditions across the boundary surface joining 
the interior space-time and the exterior space-time described by the Vaidya|4| metric 

ds% = - (l - ^M\ dv 2 - 2dvdf + f 2 d{dd 2 + sin 2 dd<p 2 ). (22) 

In Eq. fT22|) . v denotes the retarded time and m(v) represents the total mass of the collapsing 
star. The junction conditions are determined by assuming a time-like 3-surface S which 



separates the interior and the exterior manifolds 



5]. Continuity of the metric space-times 



((rfs^s = (ds 2 ^)^ = <is|) and the extrinsic curvatures (K^ = K^) across the surface E, 
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then yield the following matching conditions linking smoothly the interior (r < r%) and the 
exterior (r > r-%) space-times across the boundary: 



, 2m n dr\ 2 , 
A Q (r^)dt=-- ( 1- — + 2— J dv, 



rsf(t) =r E (f), 
rf(t) 



m[v) 



[Prk = [qf(t)B ] 
m(r,t) = m(v). 



r 2 ( rR 

R 2 + r 2+ \ A 



^j • 



(23) 
(24) 

(25) 

(26) 

(27) 



From Eq. ( 125]) and (127]) . the mass enclosed within the boundary surface at any instant t 
within a boundary surface r <r^ can be written as 



m(r,t) 



rf(t) 



r f 
IV - r' ~ [a~o 



Combining Eqs. (JBJ), © and (126"]), together within the condition (p r ) s (r 
deduce the surface equation governing the collapsing matter in the form 



(28) 



r s ) = 0, we 



2ff + f 2 - 2nf = 0, 



where, we have defined 



n 



4, 

Bo 



(29) 



(30) 



Note that for agiven initial static configuration n appears as a constant in Eq. (129]) . Fol 
lowing Bonnor|6j], we write Eq. ( 129]) as a first order differential equation 

2n 



f 



s/f 



1 



/), 



which admits a solution 



1 

t = - 

n 



Z + y7 + in(i_y7) 



(31) 



(32) 



Note that at t — > — cxd, i.e., at the onset of collapse, / = 1 and f(t) — > as t — > 0. 

We, therefore, have a complete description of the interior and exterior space-times of 
the collapsing body. In the following section, we shall analyze the impact of anisotropy by 
making use of the solutions thus obtained. 



8 



V. PHYSICAL ANALYSIS: 



To understand the nature of collapse, we assume that a star starts its collapse at t = — oo 
(i.e., f(t) = 1) and the initial static star is described by the parameters A , B and A s . 
We assume that the collapsing object has an initial radius r s (t —¥ — oo) = r s and mass 
m(r s , — oo) = m s satisfying the condition 2m s /r s < 1. If for some reason, instability develops 
inside the star, it begins to collapse. The comoving boundary surface (r/)s = r s f(t) then 
starts shrinking until it reaches its Schwarzschild horizon value [r f (£bh)] s = r sf{tbh) = 
2m(v), where t = tj,h denotes the time of formation of the black hole corresponding to the 
value of f(t) = f bh . 

From Eq. (1251) . the mass of the evolving star at any instant t within the boundary radius 
rs may be written as 



"'./ + ^(i-V7) s 



4 2 
where, mass of the initial static star has the form 

m s = Airr 2 p s dr = , - s — -- 

J Hs 2(r s 2 + # 2 ) 

Consequently, the condition (r/^)^ = 2m{y) yields 

2 
2nr 

fbh 



(33) 
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2llT 
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2m s 



(35) 



and the time of black hole formation is obtained as 

fbh 



t 



bh 



1 

n 



+ ^/f bh + \n(l- y/fi 



bh ) 



(36) 



The model parameters (namely, R, P, Q) involving the initial static star are determined 
using the following boundary conditions: 



A n (r 



OV'sJ 



B (r s ) 
(Pr)s{r s ) 



2m f 



2m, 



1/2 



-1/2 



0. 



(37) 

(38) 
(39) 



where we have matched the static interior space-time to the Schwarzschild exterior and 
imposed the condition that the radial pressure ((p r ) s ) must vanish at the boundary. 



TABLE I. Data for collapsing systems having different values of the anisotropic parameter a. All 
the configurations start collapsing with initial masses m s = 3.25 Mq and radii r s = 20 km. 

a P Q R (km) n f bh t bh r bh (km) m bh (M ) 



0.9 


0.2848 


-0.6341 


20.8427 


0.0244 


0.2298 


-2.3914 


4.596 


2.2980 


0.6 


0.3800 


-0.4994 


20.8427 


0.0163 


0.2298 


-3.5816 


4.596 


2.2980 


0.3 


0.4065 


-0.4619 


20.8427 


0.0142 


0.2298 


-4.1074 


4.596 


2.2980 





0.42067 


-0.4418 


20.8427 


0.0131 


0.2298 


-4.4578 


4.596 


2.2980 


-0.3 


0.4300 


-0.4287 


20.8427 


0.0124 


0.2298 


-4.7239 


4.596 


2.2980 


-0.6 


0.4367 


-0.4192 


20.8427 


0.0118 


0.2298 


-4.9401 


4.596 


2.2980 


-0.9 


0.4418 


-0.4120 


20.8427 


0.0114 


0.2298 


-5.1230 


4.596 


2.2980 


-1.5 


0.4493 


-0.4012 


20.8427 


0.0108 


0.2298 


-5.4231 


4.596 


2.2980 



Now, to get an insight about the effects of anisotropy, we consider different initial static 
configurations, both isotropic a = and anisotropic (a ^ 0). We consider different initially 
static stellar configurations of identical initial masses and radii (we have assumed m s = 
3.25 Mq and radius r s = 20 km) with different values of the anisotropic parameter a. 
Using numerical procedures, we have calculated the corresponding model parameters and 
also evaluated the time of formation of black holes t bh and radius (r bh = r s f bh ) and mass 
(m b h) of the black hole formed. Our results have been compiled in Table [B Variations of 
f(t) for different choices of the anisotropic parameter a have been shown in Fig. [U We have 
also calculated the rate of collapse in our model for different anisotropic parameters. The 
collapse rate, in our model, is obtained as 



3/ _ 6n(y/7- l 
;/3 A f ~ Aof^T 



e=*4 = -rr= 7;^ ■ ( 4 o) 



The collapse rate turns out to be 9 = -0.101464, - 0.102096, - 0.100749 for a = 
0, 0.9, — 0.9, respectively. Since, the collapsing object contracts in size as time evolves G, 
in our construction, turns out to be negative. However, from the absolute values of 0, we 
note that for a positive anisotropy (p r > p t ) the collapse rate increases as compared to an 
isotropic star while for a negative value of a (pt > p r ), the rate decreases. From Tabled! we 
note that 
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FIG. 1. Variation of f(t) for different anisotropic parameters. 



Case I: when a > (implying p r > p t ), the horizon is formed at a faster rate as 
compared to a = 0, i.e., isotropic case. 

Case II: When a < (implying p t > p r ), the horizon is formed at a slower rate as 
compared to a = 0, i.e., isotropic case. 



Similar observations may be found in Ref. ll|. However, we note that mass and radius of 
the collapsed configuration do not depend on anisotropy in our formulation. 



VI. THERMAL BEHAVIOUR 



To analyze the impact of anisotropy on the evolution of temperature of the collapsing 
system, we use the relativistic Maxwell-Cattaneo relation for temperature governing the 
heat transport [20, [58|, |59( given by 



r{g afi + u a u p )u b q p . s + q a = -K(g afl + u a u p )[Tp + Tu p ] 



(41) 



where k(> 0) is the thermal conductivity and r(> 0) is the relaxation time. For the line 
element (pQ), Eq. (J4T1) reduces to 



T±{ q fB )WfAoB = -K±±(A > n 



(42) 
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FIG. 2. Evolution of surface temperature for different anisotropic parameters. 



Following an earlier work 45| , we write the relativistic Fourier heat transport equation by 
setting r = in ( |4"2]) . For r = 0, combining Eqs. (jSJ) and f|42|) . we get 



8ttk(A T) 



, _ 24,/ 

" A f' 



(43) 



Let us now assume that the thermal conductivity varies as k = 7T W , where 7 and u are 
constants. Eq. (l4"3l) . then yields 



8tt 7 (A)T)' 



2A / T-< J 

An 



2n(y/7-i; 

/v7 



(44) 



where we have used Eq. (!3T|) . Integrating the above equation, we get 



/--■^ "^■"iV ^H-roft), 



27T7/V7 V A 
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(45) 



To get a simple estimate of temperature evolution, we set u = 0, 7 = 1 and T (t) = and 
evaluate the surface temperature at any instant as 



T(rx,t) 



n 



(V?-l) /In A 



2tt/v7 v A ; s 



(46) 



Time evolution of the surface temperature for different anisotropic parameter a has been 
shown in Fig. (J5J), where we have used the data from Table [0 
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VII. DISCUSSIONS 

We have developed a minimalistic and basic framework of a gravitationally collapsing 
system which has allowed us to examine the impact of pressure anisotropy explicitly. In 
our construction, we have ignored the effects of various other factors relevant for collapsing 
systems such as shear, viscosity, charge etc. Moreover, the line element describing the 
interior space-time has been assumed to be spherically symmetric where the metric functions 
have been chosen to be separable in its variables. Though formulation of a more general 
framework is always preferred to describe a realistic situation, the simple model developed 
here provides a mechanism to capture the impact of anisotropy on gravitational collapse 
successfully. Though, as pointed out in [42[, effects of factors like shear, charge etc., on 
self-bound systems can be absorbed by considering the system to be anisotropic, in general, 
we intend to formulate a more general framework so as to examine the combined impacts of 
various factor relevant for gravitationally collapsing systems. These issues will be taken up 
elsewhere. 



ACKNOWLEDGMENTS 

RS is thankful to the Inter University Centre for Astronomy and Astrophysics (IUCAA), 
Pune, India, where a part of this work was carried out under its Visiting Research Associ- 
ateship Programme. 



[1] P. S. Joshi & D. Malafarina, Int. J. Mod. Phys. D 20, 2641, 2011. 

[2] S. Chandrasekhar, Observatory 57, 373, 1934. 

[3] J. R. Oppenheimer & H. Snyder, Phys. Rev. 56, 455, 1939. 

[4] P. C. Vaidya, Proc. Indian Acad. Sci. A33, 264, 1951. 

[5] N. O. Santos, Mon. Not. R. Astron. Soc. 216, 403, 1985. 

[6] W. B. Bonnor, A. K. G. de Oliveira & N. O. Santos, Phys. Rep. 181, 269, 1989. 

[7] A. K. G. de Oliveira & N. O. Santos, Astrophys. J. 312, 640, 1987; Mon. Not. R. Astron. Soc. 
216, 1001, 1988. 

13 



[8] A. K. G. de Oliveira, C. A. Kolassis k N. O. Santos, Mon. Not. R. Astron. Soc. 231, 1011, 

1988. 
[9] L. Herrera, J. Ospino k A. di Prisco, Phys. Rev. D 77, 027502, 2008. 
[10] L. Herrera, N. O. Santos k A. Wang, Phys. Rev. D 78, 084026, 2008. 
[11] L. Herrera, A. di Prisco, J. L. Herndndez-Pastora k N. O. Santos, Phys. Lett. A 237, 113, 

1998. 
[12] L. Herrera k N. O. Santos, Mon. Not. R. Atron. Soc. 287, 161, 1997. 
[13] L. Herrera, A. di Prisco, J. Martin, J. Ospino, N. O. Santos k O. Troconis, Phys. Rev. D 69, 

084026, 2004. 
[14] L. Herrera k W. Barreto, Int J. Mod. Phys. D 20, 1265, 2011. 
[15] L. Herrera k N. O. Santos, Phys. Rep. 286, 53, 1997. 
[16] L. Herrera, J. Martin k J. Ospino, J. Math. Phys. 43,4889 2002. 
[17] L. Herrera, Int. J. Mod. Phys. D 20, 1689, 2011. 
[18] A. di Prisco, L. Herrera, G. L. Denmat, M. A. H. MacCallum k N. O. Santos, Phys. Rev. D 

76, 064017, 2007. 
[19] A. di Prisco, L. Herrera k V. Varela, Gen. Relativ. Grav. 29, 1239, 1997. 
[20] J. Martinez, Phys. Rev. D 53, 6921, 1996. 
[21] G. Pinheiro k R. Chan, Gen. Relativ. Grav. 40,2149, 2008; Gen. Relativ. Grav. 43, 1451, 

2011; Gen. Relativ. Grav. 45,243, 2013. 
[22] R. Chan, Mon. Not. R. Astro. 316, 588, 2000; Astron. & Astrophys. 368, 325, 2001. 
[23] R. Chan, M. F. A. da Silva k J. F. V. da Rocha, J. Math. Phys. D 12, 347, 2003. 
[24] R. Chan, L. Herrera & N. O. Santos Mon. Not. R. Astro. 267, 637, 1994. 
[25] R. Tikekar & L. K. Patel, Pramana-j. of physics 39, 17, 1992. 
[26] S. D. Maharaj k M. Govender, Pramana-j. of physics 54, 715, 2000. 
[27] M. Sharif k K. Iqbal, Mod. Phys. Lett. A 24, 1533, 2009. 

[28] M. Sharif & A. Siddiqe, Mod. Phys. Lett. A 25, 2831, 2010; Gen. Relativ. Grav. 43, 73, 2011. 
[29] M. Sharif k G. Abbas, J. Phys. Soc. Jpn. 80, 104002, 2011; Astrophys. Space Sci. 335, 515, 

2011. 
[30] M. Sharif & F. Sundas, Gen. Relativ. Grav. 43, 127, 2011. 
[31] M. Sharif & S. Fatima, Gen. Relativ. Grav. 43, 127, 2011. 
[32] W. Barreto, B. Rodriguez, L. Rosales k O. Serrano, Gen. Relativ. Grav. 39, 23, 2007. 

14 



[33] C. Ghezzi, Phys. Rev. D 72, 104017, 2005. 

[34] R. Goswami k P. S. Joshi, Class. Quantum Grav. 21, 3645, 2004; Class. Quantum Grav. 19, 

0129, 2002; Phys. Rev D 69, 027502, 2004. 
[35] S. Chakraborty k T. Bandyopadhyay, Int. J. Mod. Phys. D 17, 1271, 2008. 
[36] D. M. Eardley k L. Smarr, Phys. Rev. D 19, 2239, 1979. 
[37] F. C. Mena, B. C. Nolan k R. Tavakol, Phys. Rev. D 70, 084030, 2004. 
[38] S. Thirukkanesh k S. D. Maharaj, Math. Meth. Appl. Sci. 32, 684, 2009. 
[39] K. S. Govinder, M. Govender k R. Maartens, Mon. Not. R. Astro. Soc. 299, 809, 1998. 
[40] D. Schafer k H. F. Goenner, Gen. Relativ. Grav. 42, 2119, 2000. 
[41] B. V. Ivanov, Gen Relativ. Grav. 44, 1835, 2012; Int. J. Mod. Phys. D 20, 319, 2011; Int. J. 

Mod. Phys. A 25, 3975, 2010. 
[42] B. V. Ivanov, Int. J. Theor. Phys. 49, 1236, 2010. 
[43] B. V. Ivanov, Phys. Rev. D 65, 104011, 2002. 
[44] S. Sarwe k R. Tikekar, Int. J. Mod. Phys. D 19, 1889 (2010). 
[45] R. Sharma k R. Tikekar, Gen. Relativ. Grav. 44, 2503, 2012. 
[46] R. Sharma k R. Tikekar, Pramana- j. of physics 79, 501, 2012. 
[47] R. L. Bowers and E. P. T. Liang, Astrophys. J. 188, 657, 1974. 
[48] M. K. Mak k T. Harko, Proc. Roy. Soc. Lond. A 459, 393, 2003. 
[49] T. Harko k M. K. Mak, Annalen Phys. 11, 3, 2002. 
[50] K. Dev k M. Gleiser, Int. J. Mod. Phys. D 13, 1389, 2004. 
[51] S. Karmakar, S. Mukherjee, R. Sharma and S. D. Maharaj, Pramana-.!. of Phys. 68, 881, 

2007. 
[52] M. R. Finch k J. E. F. Skea, Class. Quantum. Grav. 6, 467, 1989. 
[53] S. Hansraj k S. D. Maharaj, Int. J. Mod. Phys. D 15, 1311, 2006. 
[54] P. C. Vaidya k R. Tikekar, J. Astrophys. Astron. 3, 325, 1982. 
[55] R. Tikekar k K. Jotania, Gravitation and cosmology 15, 129, 2009. 
[56] R. Tikekar k K. Jotania, Pramana- J. of Phys. 68, 397, 2007. 
[57] K. Jotania k R. Tikekar, Int. J. Mod. Phys. D 15, 1175, 2006. 
[58] W. Israel k J. Stewart, Ann. Phys. 118, 341, 1979. 
[59] R. Maartens, Class. Quantum Grav. 12, 1455, 1995. 



15 



